Diversity of bacterial type II toxin–antitoxin systems: a comprehensive search and functional analysis of novel families

Type II toxin–antitoxin (TA) systems are generally composed of two genes organized in an operon, encoding a labile antitoxin and a stable toxin. They were first discovered on plasmids where they contribute to plasmid stability by a phenomenon denoted as ‘addiction’, and subsequently in bacterial chromosomes. To discover novel families of antitoxins and toxins, we developed a bioinformatics approach based on the ‘guilt by association’ principle. Extensive experimental validation in Escherichia coli of predicted antitoxins and toxins increased significantly the number of validated systems and defined novel toxin and antitoxin families. Our data suggest that toxin families as well as antitoxin families originate from distinct ancestors that were assembled multiple times during evolution. Toxin and antitoxin families found on plasmids tend to be promiscuous and widespread, indicating that TA systems move through horizontal gene transfer. We propose that due to their addictive properties, TA systems are likely to be maintained in chromosomes even though they do not necessarily confer an advantage to their bacterial hosts. Therefore, addiction might play a major role in the evolutionary success of TA systems both on mobile genetic elements and in bacterial chromosomes

Ongoing diphtheria outbreak in Yemen: a cross-sectional and genomic epidemiology study.

BACKGROUND: An outbreak of diphtheria, declared in Yemen in October, 2017, is ongoing. We did a cross-sectional study to investigate the epidemiological, clinical, and microbiological features of the outbreak. METHODS: Probable cases of diphtheria that were defined clinically and recorded through a weekly electronic diseases early warning system (from 2017, week 22, to 2020, week 17) were used to identify trends of the outbreak (we divided the epidemic into three time periods: May 29, 2017, to June 10, 2018; June 11, 2018, to June 3, 2019; and June 4, 2019, to April 26, 2020). We used the line list of diphtheria reports for governorate-level descriptions. Vaccination coverage was estimated using the 2017 and 2018 annual reports by the national Expanded Programme on Immunization. To confirm cases biologically, Corynebacterium diphtheriae was isolated and identified from throat swabs using standard microbiological culture and identification procedures. We assessed differences in the temporal and geographical distributions of cases, including between different age groups. For in-depth microbiological analysis, tox gene and species-specific rpoB real-time PCR, Illumina genomic sequencing, antimicrobial susceptibility analysis (disk diffusion, E-test), and the Elek diphtheria toxin production test were done on confirmed cases. We used genomic data for phylogenetic analyses and to estimate the nucleotide substitution rate. FINDINGS: The Yemen diphtheria outbreak affected almost all governorates (provinces), with 5701 probable cases and 330 deaths recorded up to April 26, 2020. We collected clinical data for 888 probable cases with throat swab samples referred for biological confirmation, and genomic data for 42 positive cases, corresponding to 43 isolates (two isolates from one culture were included due to distinct colony morphologies). The median age of patients was 12 years (range 0·2-80). The proportion of cases in children aged 0-4 years was reduced during the second time period, after a vaccination campaign, compared with the first period (19% [95% CI 18-21] in the first period vs 14% [12-15] in the second period, p<0·0001). Among 43 tested isolates, 39 (91%) produced the diphtheria toxin and two had low level (0·25 mg/L) antimicrobial resistance to penicillin. We identified six C diphtheriae phylogenetic sublineages, four of which are genetically related to isolates from Saudi Arabia, Eritrea, and Somalia. Inter-sublineage genomic variations in genes associated with antimicrobial resistance, iron acquisition, and adhesion were observed. The predominant sublineage (30 [70%] of 43 isolates) was resistant to trimethoprim and was associated with unique genomic features, more frequent neck swelling (p=0·0029) and a younger age of patients (p=0·060) compared with the other sublineages. Its evolutionary rate was estimated at 1·67 × 10-6 substitutions per site per year, placing its most recent common ancestor in 2015, and indicating silent circulation of C diphtheriae in Yemen before the outbreak was declared. INTERPRETATION: In the Yemen outbreak, C diphtheriae shows high phylogenetic, genomic, and phenotypic variation. Laboratory capacity and real-time microbiological monitoring of diphtheria outbreaks need to be scaled up to inform case management and transmission control of diphtheria. Catch-up vaccination might have provided some protection to the targeted population (children aged 0-4 years). FUNDING: National Centre of the Public Health Laboratories (Yemen), Institut Pasteur, and the French Government Investissement d'Avenir Programme. TRANSLATION: For the Arabic translation of the abstract see Supplementary Materials section

Key components of the eight classes of type IV secretion systems involved in bacterial conjugation or protein secretion

Conjugation of DNA through a type IV secretion system (T4SS) drives horizontal gene transfer. Yet little is known on the diversity of these nanomachines. We previously found that T4SS can be divided in eight classes based on the phylogeny of the only ubiquitous protein of T4SS (VirB4). Here, we use an ab initio approach to identify protein families systematically and specifically associated with VirB4 in each class. We built profiles for these proteins and used them to scan 2262 genomes for the presence of T4SS. Our analysis led to the identification of thousands of occurrences of 116 protein families for a total of 1623 T4SS. Importantly, we could identify almost always in our profiles the essential genes of well-studied T4SS. This allowed us to build a database with the largest number of T4SS described to date. Using profile-profile alignments, we reveal many new cases of homology between components of distant classes of T4SS. We mapped these similarities on the T4SS phylogenetic tree and thus obtained the patterns of acquisition and loss of these protein families in the history of T4SS. The identification of the key VirB4-associated proteins paves the way toward experimental analysis of poorly characterized T4SS classes

Phylogenomics of Unusual Histone H2A Variants in Bdelloid Rotifers

Rotifers of Class Bdelloidea are remarkable in having evolved for millions of years, apparently without males and meiosis. In addition, they are unusually resistant to desiccation and ionizing radiation and are able to repair hundreds of radiation-induced DNA double-strand breaks per genome with little effect on viability or reproduction. Because specific histone H2A variants are involved in DSB repair and certain meiotic processes in other eukaryotes, we investigated the histone H2A genes and proteins of two bdelloid species. Genomic libraries were built and probed to identify histone H2A genes in Adineta vaga and Philodina roseola, species representing two different bdelloid families. The expressed H2A proteins were visualized on SDS-PAGE gels and identified by tandem mass spectrometry. We find that neither the core histone H2A, present in nearly all other eukaryotes, nor the H2AX variant, a ubiquitous component of the eukaryotic DSB repair machinery, are present in bdelloid rotifers. Instead, they are replaced by unusual histone H2A variants of higher mass. In contrast, a species of rotifer belonging to the facultatively sexual, desiccation- and radiation-intolerant sister class of bdelloid rotifers, the monogononts, contains a canonical core histone H2A and appears to lack the bdelloid H2A variant genes. Applying phylogenetic tools, we demonstrate that the bdelloid-specific H2A variants arose as distinct lineages from canonical H2A separate from those leading to the H2AX and H2AZ variants. The replacement of core H2A and H2AX in bdelloid rotifers by previously uncharacterized H2A variants with extended carboxy-terminal tails is further evidence for evolutionary diversity within this class of histone H2A genes and may represent adaptation to unusual features specific to bdelloid rotifers

The repertoire of ICE in prokaryotes underscores the unity, diversity, and ubiquity of conjugation

Horizontal gene transfer shapes the genomes of prokaryotes by allowing rapid acquisition of novel adaptive functions. Conjugation allows the broadest range and the highest gene transfer input per transfer event. While conjugative plasmids have been studied for decades, the number and diversity of integrative conjugative elements (ICE) in prokaryotes remained unknown. We defined a large set of protein profiles of the conjugation machinery to scan over 1,000 genomes of prokaryotes. We found 682 putative conjugative systems among all major phylogenetic clades and showed that ICEs are the most abundant conjugative elements in prokaryotes. Nearly half of the genomes contain a type IV secretion system (T4SS), with larger genomes encoding more conjugative systems. Surprisingly, almost half of the chromosomal T4SS lack co-localized relaxases and, consequently, might be devoted to protein transport instead of conjugation. This class of elements is preponderant among small genomes, is less commonly associated with integrases, and is rarer in plasmids. ICEs and conjugative plasmids in proteobacteria have different preferences for each type of T4SS, but all types exist in both chromosomes and plasmids. Mobilizable elements outnumber self-conjugative elements in both ICEs and plasmids, which suggests an extensive use of T4SS in trans. Our evolutionary analysis indicates that switch of plasmids to and from ICEs were frequent and that extant elements began to differentiate only relatively recently. According to the present results, ICEs are the most abundant conjugative elements in practically all prokaryotic clades and might be far more frequently domesticated into non-conjugative protein transport systems than previously thought. While conjugative plasmids and ICEs have different means of genomic stabilization, their mechanisms of mobility by conjugation show strikingly conserved patterns, arguing for a unitary view of conjugation in shaping the genomes of prokaryotes by horizontal gene transfer

Origine et évolution des systèmes toxine-antitoxine de classe II

Les systèmes toxine-antitoxine (TA) sont composés de deux gènes organisés en opéron retrouvés chez la quasi-totalité des bactéries ainsi que chez les archées. Ils ont été découverts sur des plasmides, où ils induisent une tuerie post-segregationnelle (PSK). En effet, l’antitoxine est instable car elle est dégradée par une protéase ATP dépendante. Lors de la réplication, si un plasmide portant un système TA n’est pas transmis à la cellule fille, celle-ci recevra tout de même une partie du cytoplasme de la cellule mère qui contenait des protéines de toxine et d’antitoxine. Cette dernière étant instable, et sans synthèse de novo, la toxine va se retrouver libre d’affecter sa cible cellulaire. Un système TA crée donc une addiction au plasmide qui le porte, stabilisant celui-ci.Les systèmes TA se retrouvent également, dans un nombre de copies variable, au sein du chromosome. Dans ce cas, plusieurs hypothèses existent quant à leur fonction, comme la mort cellulaire programmée, la réponse au stress, la stabilisation de régions génomiques non essentielles, ou l’anti-addiction au plasmide.L’origine et l’évolution des systèmes TA restent inconnues, alors qu’ils présentent des aspects intrigants. En effet, les toxines CcdB et MazF ont des séquences et des activités toxiques très différentes, malgré une structure proche ;les toxines RelE et ParE présentent des séquences proches, mais des fonctions différentes. À l’heure actuelle, les deux hypothèses concernant l’origine évolutive des systèmes TA sont soit qu’ils seraient tous issus d’un ancêtre commun, soit qu’ils auraient été réinventés plusieurs fois au cours de l’évolution, à partir d’un nombre limité de gènes.Au cours de cette thèse, nous avons abordé la question de l’origine et l’évolution des systèmes TA de plusieurs manières :par la reconstruction de phylogénies et de séquences ancestrales, et par l’étude d’un système chromosomique particulier au sein de nombreuses souches d’Escherichia coli. Enfin, nous nous avons décidé d’analyser le contexte génomique des systèmes TA chromosomiques. Ces travaux ont notamment permis de mieux comprendre l’évolution des systèmes TA, et de conforter l’hypothèse selon laquelle ils seraient des éléments égoïstes, évoluant au sein des génomes sans gain d’aptitude pour l’hôte.Doctorat en Sciencesinfo:eu-repo/semantics/nonPublishe

Evolutionary history of hepatitis C virus genotype 5a in France, a multicenter ANRS study

The epidemic history of HCV genotype 5a is poorly documented in France, where its prevalence is very low, except in a small central area, where it accounts for 14.2% of chronic hepatitis C cases. A Bayesian coalescent phylogenetic investigation based on the E1 envelope gene and a non-structural genomic segment (NS3/4) was carried out to trace the origin of this epidemic using a large sample of genotype 5a isolates collected throughout France. The dates of documented transmissions by blood transfusion were used to calibrate five nodes in the phylogeny. The results of the E1 gene analysis showed that the best-fitting population dynamic model was the expansion growth model under a relaxed molecular clock. The rate of nucleotide substitutions and time to the most recent common ancestors (tMRCA) of genotype 5a isolates were estimated. The divergence of all the French HCV genotype 5a strains included in this study was dated to 1939 [95% HPD: 1921-1956], and the tMRCA of isolates from central France was dated to 1954 [1942-1967], which is in agreement with epidemiological data. NS3/4 analysis provided similar estimates with strongly overlapping HPD values. Phylodynamic analyses give a plausible reconstruction of the evolutionary history of HCV genotype 5a in France, suggesting the concomitant roles of transfusion, iatrogenic route and intra-familial transmission in viral diffusion. © 2011 Elsevier B.V.0SCOPUS: ar.jSCOPUS: ar.jinfo:eu-repo/semantics/publishe

Simulation data for the estimation of numerical constants for approximating pairwise evolutionary distances between amino acid sequences

International audienceEstimating the number of substitution events per site that have occurred during the evolution of a pair of amino acid sequences is a common task in phylogenetics and comparative genomics that often requires quite slow maximum-likelihood procedures when taking into account explicit evolutionary models. Data presented in this article are large sets of numbers of substitution events and associated numbers of observed differences between pairs of aligned amino acid sequences that have been generated through a simulation procedure of sequence evolution under a broad range of evolutionary models. These data are available at https://zenodo.org/record/2653704 (doi:10.5281/zenodo.2653704). They are accompanied in this paper by figures showing the strong relationship between the corresponding evolutionary and uncorrected distances, as well as estimated numerical constants that determine non-linear functions that fit the simulated data. These numerical constants can be useful to quickly estimate pairwise evolutionary distances directly from uncorrected distances between aligned amino acid sequences

Automated discovery and phylogenetic analysis of new toxin-antitoxin systems.

Although often viewed as elements "at the service of" bacteria, plasmids exhibit replication and maintenance mechanisms that make them purely "selfish DNA" candidates. Toxin-antitoxin (TA) systems are a spectacular example of such mechanisms: a gene coding for a cytotoxic stable protein is preceded by a gene coding for an unstable antitoxin. The toxin being more stable than the antitoxin, absence of the operon causes a reduction of the amount of the latter relative to the amount of the former. Thus, a cell exhibiting a TA system on a plasmid is 'condemned' either not to loose it or to die.Journal ArticleResearch Support, Non-U.S. Gov'tSCOPUS: ar.jinfo:eu-repo/semantics/publishe